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Abstract 

Neutral or charged polymers that are densely end-grafted to surfaces form brush- 
like structures and are highly stretched under good-solvent conditions. We discuss and 
compare relevant results from scaling models, self-consistent-field methods and MD 
simulation techniques and concentrate on the conceptual simple case of planar sub- 
strates. For neutral polymers the main quantity of interest is the brush height and the 
polymer density profile, which can be well predicted from self-consistent calculations 
and simulations. Charged polymers (polyelectrolytes) are of practical importance since 
they are soluble in water. Counterion degrees of freedom determine the brush behavior 
in a decisive way and lead to a strong and non-linear swelling of the brush. 
Keywords: Brushes, polyelectrolytes, scaling theory, self-consistent field theory, sim- 
ulation techniques 
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Legend of Symbols 

a: Kuhn length or effective monomer size 

d: height of counterion layer 

/: fractional charge of the chain < / < 1 

F: free energy in units of ksT (per chain or unit area) 

h: height of brush 

fc^T: thermal energy 

L: contour length of a chain 

£b = / {AneksT): Bjerrum length 

N: polymerization index 

R: end-to-end polymer chain radius 

Rq: end-to-end radius of an ideal polymer 

Rp: Flory radius of a self-avoiding chain 

V2: 2nd virial coefficient of monomers in solution 

K~^: Debye-Hiickel screening length 

u: Flory exponent for the polymer size 

11: osmotic pressure in units of ksT 

Pa- grafting density of a polymer brush 

p{z): monomer density at distance z from grafting surface 

(7: Lennard-Jones diameter in simulation 

1 Introduction 

Polymers are long, chain-like molecules that consist of repeating subunits, the so-called 
monomers[l, 2, 3, 4]. In many situations, all monomers of a polymer are alike, showing 
for example the same tendency to adsorb to a substrate[5]. For industrial applications, 
one is often interested in end-functionalized polymers that are attached with one end only 
to the substrate[6, 7]. Industrial interest comes from the need to stabilize particles and 
surfaces against fiocculation. In end-grafted polymer structures the stabilization power is 
greatly enhanced as compared with adsorbed layers of polymers, where each monomer is 
equally attracted to the substrate. The main reason is that bridging of polymers between 
two approaching surfaces and creation of polymer loops on the same surface is very frequent 
in the case of polymer adsorption and eventually leads to attraction between two particle 
surfaces and thus destabilization. This does not occur if the polymer is grafted by its end 
to the surface and the monomers are chosen such that they do not particularly adsorb to 
the surface. In biology, brush-like polymer structures are encountered as the surface coating 
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of endothelial cells and regulate the adsorption and migration of various particles and bio- 
molecules from the blood stream to the vascular cells. 

Experimentally, two basic ways of building a grafted polymer layer can be distinguished: 
In the first, the polymerization is started from the surface with some suitably chosen surface- 
linked initiator. The advantage of this grafting- from procedure is that only monomers have to 
diffuse through the forming brush layer and thus the reaction kinetics is fast. In the second 
route one attaches polymers with special end-groups that act as anchors on the surface. 
This grafting-to procedure is subject to slow kinetics during the formation stage since whole 
polymers have to diffuse through the natant grafting layer, but benefits from a somewhat 
better control over the brush constitution and chemical composition. One distinguishes 
physical adsorption of end-groups that favor the substrate, for example zwitter-ionic end- 
groups attached to poly-styrene chains that lead to binding to mica in organic solvents 
such as toluene or xylene [8]. A stronger and thus more stable attachment is possible with 
cowa/en^/?/ end- grafted chains, for example poly-dimethylsiloxane chains which carry hydroxyl 
end groups and undergo condensation reactions with silanols of a silica surface [9]. One can 
also employ a suitably chosen diblock copolymer where one block adsorbs on the substrate 
while the other is repelled from it [10]. An example is furnished by polystyrene-poly (vinyl- 
pyridine) (PS-PVP) diblocks in the selective solvent toluene, which is a bad solvent for the 
PVP block and promotes strong adsorption on a quartz substrate, but acts as a good solvent 
for the PS block and disfavors its adsorption on the substrate [11]. In a slight modification 
one uses diblock copolymers that are anchored at the liquid-air [12, 13] or at a liquid-liquid 
interface of two immiscible liquids [14]. This scenario offers the advantage that the grafting 
density can (for the case of strongly anchored polymers) be varied by lateral compression (like 
a Langmuir mono-layer) and that the lateral surface pressure can be directly measured. The 
lateral pressure is an important thermodynamic quantity and allows detailed comparison 
with theoretical predictions. A well studied example is that of a diblock copolymer of 
polystyrene - polyethylene oxide (PS-PE0)[13]. The PS block is shorter and functions as 
an anchor at the air/water interface because it is immiscible in water. The PEO block is 
miscible in water but because of attractive interaction with the air/water interface it forms 
a quasi-two dimensional layer at very low surface coverage. As the surface pressure increases 
and the area per polymer decreases, the PEO block is expelled from the surface and forms 
a polymer brush. 

In this theoretical chapter we simplify the discussion by assuming that the polymers are 
irreversibly grafted at one of their chain ends to the substrate, we only mention in passing 
a few papers on the kinetics of grafting. The substrate is assumed to be solid, planar and 
impenetrable to the the polymer monomers, we briefiy cite some results obtained for curved 
substrates. We limit the discussion to good solvent conditions and neglect any attractive 
interactions between the polymer chains and the surface. Charged polymers are interesting 
from the application point of view, since they allow for water-based formulations of organic 
substances which are advantageous for economical and ecological reasons. Recent years 
have seen a tremendous research activity on charged polymers in bulk[15, 16, 17, 18] and at 
interfaces[5, 19]. We therefore treat neutral brushes as well as charged ones. 

The characteristic parameter for brush systems is the anchoring or grafting density pa, 
which is the inverse of the average area available for each polymer at the surface. For 
small grafting densities, pa < Pa, the polymer chains will be far apart from each other 
and hardly interact, as schematically shown in Fig. la. The polymers in this case form well 
separated mushrooms at the surface. The grafting density at which chains just start to 
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overlap is determined by p* ~ where R is the typical radius or size of a chain. In good 
solvent conditions (that is for swollen chains), the chain radius follows the Flory prediction 
Rp ^ aN^^^ where N is the polymerization index or monomer number of the chain, and a is 
a characteristic microscopic length scale of the polymer that incorporates monomer size as 
well as backbone stiffness. The crossover grafting density for a polymer under good-solvent 
conditions follows as p* ~ a~'^N~^/^ . For large grafting densities, p > p*, the chains are 
strongly overlapping. This situation is depicted in Fig. lb. Since we assume the solvent to 
be good, monomers repel each other. The lateral separation between the polymer coils is 
fixed by the grafting density, so that the polymers extend away from the grafting surface in 
order to avoid each other. The resulting structure is called a polymer brush, with a vertical 
height h which greatly exceeds the unperturbed coil radius R [20, 21, 22]. Similar stretched 
structures occur in many other situations, such as diblock copolymer melts in the strong 
segregation regime [23, 6], or star polymers under good solvent conditions [24]. Theory is 
mostly concerned with predicting the layer height h, but also the detailed monomer density 
profile and resulting forces of lateral or vertical brush compression as a function of the various 
system parameters. 

The understanding of grafted polymer systems progressed substantially with the advent of 
experimental techniques such as: surface-force-balance [8], small-angle-neutron scattering [9] , 
neutron[ll, 25] and X-ray[26] diffraction, and ellipsometry[14]. Of equal merit was the 
advancement in the theoretical methodology ranging from field theoretical methods and 
scaling arguments to numerical simulations, which will be amply reviewed in this chapter. 

2 Polymer basics 

The main parameters used to describe a polymer chain are the polymerization index N, 
which counts the number of repeat units or monomers along the chain, and the size of one 
monomer or the distance between two neighboring monomers. The monomer size ranges from 
a few Angstroms for synthetic polymers to a few nanometers for biopolymers. The simplest 
theoretical description of fiexible chain conformations is achieved with the so-called freely- 
jointed chain (FJC) model, where a polymer consisting oi N + 1 monomers is represented 
by A'^ bonds defined by bond vectors with j — 1,. . . N. Each bond vector has a fixed 
length \rj\ = a corresponding to the Kuhn length, but otherwise is allowed to rotate freely 
and independently of its neighbors, as is schematically shown in Fig. 2 (top). This model 
of course only gives a coarse-grained description of real polymer chains, but we will later 
see that by a careful interpretation of the Kuhn length a and the monomer number N, an 
accurate description of the large-scale properties of real polymer chains is possible. The main 
advantage is that due to the simplicity of the FJC model, many interesting obscrvables (such 
as chain size or distribution functions) can be calculated with relative ease. We demonstrate 
this by calculating the mean end-to-end radius of such a FJC polymer. Fixing one of the 
chain ends at the origin, the position of the {k + l)-th monomer is given by the vectorial 
sum 

k 

Rfc = Er. • (1) 

Because two arbitrary bond vectors are uncorrclatcd in this simple model, the thermal 
average over the scalar product of two different bond vectors vanishes, (r^- ■ r^) = for j ^ k, 
while the mean squared bond vector length is simply given by (r^) = a^. It follows that the 
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mean squared end-to-end radius R% is proportional to the number of monomers, 

Rl = {R%) = Na"" = La, (2) 

where the contour length of the chain is given by L = Na. Rq denotes the mean end-to-end 
radius of an ideal chain, and according to Eq.(2), it scales as Rq — aN^^"^. Experimentally one 
often knows, via knowledge of the chemical structure and polymer mass, the length L and, 
via light-scattering, the radius i? of a chain. Using the above scaling results, valid for a non- 
interacting chain, the Kuhn length follows as a = R^/L and the effective monomer number 
as = L/a, which allows treatment of real chains with a complicated local conformational 
structure within the FJC model. Note that the so-determined Kuhn length a is usually 
larger than the actual monomer-size, since it takes back-bone stiffness effects into account. 
Likewise, the determined effective monomer number N is typically smaller than the actual 
(chemical) number of monomers in the chain. 

In many theoretical calculations aimed at elucidating large-scale properties, the simplifi- 
cation is carried even a step further and a continuous model is used, as schematically shown 
in Fig. 2 (bottom). In such models the polymer backbone is replaced by a continuous line 
and all microscopic details are neglected. The chain is then only characterized by its length 
L and radius R. 

2.1 Polymer swelling and collapse 

The models discussed so far describe ideal chains and do not account for interactions between 
monomers which typically consist of some short-ranged repulsion and long-ranged attraction. 
Including these interactions will give a different scaling behavior for long polymer chains. 
The end-to-end radius, R — ^ (-Rat), can be written for N ^ 1 as 

R ~ aN'' (3) 

which defines the so-called swelling exponent v. As we have seen, for an ideal polymer chain 
(no interactions between monomers), Eq.(2) implies v — 1/2. This situation is realized for 
experimental polymers at a certain temperature or solvent conditions when the attraction 
between monomers exactly cancels the stcric repulsion (which is due to the fact that the 
monomers cannot penetrate each other). This situation can be achieved in the condition 
of theta solvents. In good solvents, on the other hand, the monomer-solvent interaction is 
more attractive than the monomer-monomer interaction, in other words, the monomers try 
to avoid each other in solution. As a consequence, single polymer chains in good solvents 
have swollen spatial configurations dominated by the steric repulsion, characterized by an 
exponent u ~ 3/5, leading to the Flory radius Rp ^ aN^^^. This spatial size of a polymer 
coil is much smaller than the extended contour length L — aN but larger than the size of an 
ideal chain Rq — aN^^^, i.e. Rq <^ Rp <^ Lin the limit oo. The reason for this behavior 
is conformational entropy (which prevents full stretching of the chain) combined with the 
favorable interaction between monomers and solvent molecules in good solvents (which leads 
to a more open structure than an ideal chain). In the opposite case of bad or poor solvent 
conditions, the effective interaction between monomers is attractive, leading to collapse of 
the chains and to their precipitation from solution. In this case, the polymer volume, like any 
space filling object embedded in three-dimensional space, scales proportional to its weight 
as R^ ~ N, yielding an exponent u = 1/3. 



5 



The standard way of taking into account interactions between monomers is the Flory the- 
ory, which treats these interactions on an approximate mean- field level [1, 2, 3]. Let us first 
consider the case of repulsive interactions between monomers, which can be described by a 
positive second-virial coefficient V2 and corresponds to the good-solvent condition. For pure 
hard-core interactions and with no additional attractions between monomers, the second 
virial coefficient (which corresponds to the excluded volume) is of the order of the monomer 
volume, i.e. V2 a^- The repulsive interaction between monomers, which tends to swell the 
chain, is counteracted and balanced by the ideal chain elasticity, which is brought about by 
the entropy loss associated with stretching the chain. The origin is that the number of poly- 
mer configurations having an end-to-end radius of the order of the unperturbed end-to-end 
radius is large. These configurations are entropically favored over configurations character- 
ized by a large end-to-end radius, for which the number of possible polymer conformations 
is drastically reduced. The standard Flory theory for a fiexible chain of radius R is based on 
writing the free energy F of a chain (in units of the thermal energy ksT) as a sum of two 
terms (omitting numerical prefactors) 



where the first term is the entropic elastic energy associated with swelling a polymer chain 
to a radius R, proportional to the effective spring constant of an ideal chain, ksT/R^, and 
the second term is the second-virial repulsive energy proportional to the coefficient V2, and 
the segment density squared. It is integrated over the volume R^. The optimal radius R is 
calculated by minimizing this free energy and gives the swollen radius 

Rf - a(v2/a^)'/^iV" , (5) 

with u = 3/5. For purely steric interactions with V2 — we obtain R ~ aN'^. For 
weakly interacting monomers, < ^2 < a^, one finds that the swollen radius Eq.(5) is only 
reahzed above a minimal monomer number N ~ (^2/0^)"^ below which the chain statistics 
is unperturbed by the interaction and the scaling of the chain radius is Gaussian and given 
by Eq.(2). 

In the opposite limit of a negative second virial coefficient, V2 < 0, corresponding to the 
bad or poor solvent regime, the polymer coil will be collapsed due to attraction between 
monomers. In this case, the attraction term in the free energy is balanced by the third- virial 
term in a low-density expansion (where we assume that V3 > 0), 

F:...R'(^)\.,Ri^)\ (6) 



Minimizing this free energy with respect to the chain radius one obtains 

R^{v,/\v2\y/'N^ (7) 

with z/ = 1/3. This indicates the formation of a compact globule, since the monomer density 
inside the globule, p ~ N/R^, is independent of the chain length. The minimal chain length 
to observe a collapse behavior is ~ {v3/(i^V2y. For not too long chains and a second virial 
coefficient not too much differing from zero, the interaction is irrelevant and one obtains 
effective Gaussian or ideal behavior. It should be noted, however, that even small deviations 
from the exact theta conditions (defined by strictly V2 — 0) will lead to chain collapse or 
swelling for very long chains. 
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2.2 Charged polymers 

A polyelectrolyte (PE) is a polymer with a fraction / of charged monomers. When this 
fraction is small, / <C 1, the PE is weakly charged, whereas when / is close to unity, 
the polyelectrolyte is strongly charged. There are two common ways to control / [17]. 
One way is to polymerize a heteropolymer mixing strongly acidic and neutral monomers as 
building blocks. Upon contact with water, the acidic groups dissociate into positively charged 
protons (H"*") that bind immediately to water molecules, and negatively charged monomers. 
Although this process effectively charges the polymer molecules, the counter-ions make the 
PE solution electro-neutral on larger length scales. The charge distribution along the chain is 
quenched during the polymerization stage, and it is characterized by the fraction of charged 
monomers on the chain, /. In the second way, the PE is a weak polyacid or polybase. The 
effective charge of each monomer is controlled by the pH and the salt concentration of the 
solution[27]. Moreover, this annealed fraction depends on the local electric potential which is 
in particular important for adsorption or binding processes since the local electric potential 
close to a strongly charged surface or a second charged polymer can be very different from 
its value in the bulk solution and therefore modify the polyelectrolyte charge [28]. 

Counterions are attracted to the charged polymers via long-ranged Coulomb interac- 
tions; this physical association typically leads to a rather loosely bound counter-ion cloud 
around the PE chain. Because of this background of a polarizable and diffusive counterion 
cloud, there is a strong influence of the counterion distribution on the PE structure and 
vice versa. Counterions contribute significantly towards bulk properties, such as the osmotic 
pressure, and their translational entropy is responsible for the generally good water solubility 
of charged polymers. In addition, the statistics of PE chain conformations is governed by 
intra-chain Coulombic repulsion between charged monomers; this results in a more extended 
and swollen conformation of PE's as compared to neutral polymers and gives rise to the 
characteristically high viscosity of polyelectrolyte solutions (hence their use as viscosifiers in 
food industry). 

For polyelectrolytes, electrostatic interactions provide the driving force for their salient 
features and have to be included in any theoretical description. The reduced electrostatic 
interaction between two point- like charges can be written as qiq2v{r) where 

v{r) = Ib/t (8) 

is the Coulomb interaction between two elementary charges in units of /c^T and qi and q2 
are the valencies (or the reduced charges in units of the elementary charge e). The Bjerrum 
length £b is defined as 

^- = sir- 

where e is the medium dielectric constant. It denotes the distance at which the Coulombic 
interaction between two unit charges in a dielectric medium is equal to thermal energy (/csT). 
It is a measure of the distance below which the Coulomb energy is strong enough to compete 
with the thermal fluctuations; in water at room temperatures, one finds £b ~ 0.7 nm. 

In biological systems and most industrial applications, the aqueous solution contains 
in addition to the counterions mobile salt ions. Salt ions of opposite charge are drawn to 
the charged object and modify the counterion cloud around it. They effectively reduce or 
screen the charge of the object. The effective (screened) electrostatic interaction between 
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two charges qi and q2 in the presence of salt ions and a polarizable solvent can be written 
as qiq2VDn{r), with the Debye-Hiickel (DH) potential v-oiiiT) given (in units of ksT) by 



t'DH(r) 



= — e 
r 



(10) 



The exponential decay is characterized by the screening length k~^, which is related to the 
salt concentration Csait by /t^ = Snq'^isCsait , where q denotes the valency oi q : q salt. 
At physiological conditions the salt concentration is Cgait ~ O.IM and for monovalent ions 
{q — 1) this leads to k,~^ ^ 1 nm. This means that although the Coulombic interactions are 
long-ranged, in physiological conditions they are highly screened above length scales of a few 
nanometers, which results from multi-body correlations between ions in a salt solution. 

For charged polymers, the effective bending stiffness and thus the Kuhn length is in- 
creased due to electrostatic repulsion between monomers[29, 30, 31, 32, 33, 34, 35]. This 
effect modifies considerably not only the PE behavior in solution but also their adsorption 
characteristics [36] . 

A peculiar phenomenon occurs for highly charged PE's and is known as the Manning 
condensation of counter-ions [37, 38, 39]. This phenomenon constitutes a true phase transi- 
tion in the absence of added salt ions. For a single rigid PE chain represented by an infinitely 
long and straight cylinder with a linear charge density r larger than the threshold 



where q is the counter-ion valency, it was shown that counter-ions condense on the oppositely 
charged cylinder even in the limit of infinite system size. For a solution of stiff charged poly- 
mers this corresponds to the limit where the inter-chain distance tends to infinity. This effect 
is not captured by the linear Debye-Hiickel theory. A simple heuristic way to incorporate 
the non-linear Manning condensation is to replace the bare linear charge density r by the 
renormalized one, Trenorm — ^/{q^b), whenever isTq > 1 holds. This procedure, however, 
is not totally satisfactory at high salt concentrations [40, 41]. Also, real polymers have a 
finite length, and are neither completely straight nor in the infinite dilution limit [42]. Still, 
Manning condensation has an experimental significance for polymer solutions[43] because 
thermodynamic quantities, such as counter- ion activities [44] and osmotic coefficients [45], 
show a pronounced signature of Manning condensation. Locally, polymer segments can be 
considered as straight over length scales comparable to the persistence length. The Manning 
condition Eq. (11) usually denotes a region where the binding of counter-ions to charged 
chain sections begins to deplete the solution from free counter-ions. 

3 Neutral grafted polymers 
3.1 Scaling approach 

The scaling behavior of the brush height h can be analyzed using a Flory-like mean-field 
theory, which is a simplified version of the original Alexander theory [21] for polymer brushes. 
The stretching of the chain leads to an entropic free energy loss of /i^/ per chain, and the 
repulsive energy density due to unfavorable monomer-monomer contacts is proportional to 
the squared monomer density times the excluded- volume parameter V2. The derivation is 
thus analogous to the calculation of the Flory radius of a chain in good solvent shown in 



isrq = 1 , 



(11) 
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Section 2.1, except that now the grafting density pa plays a decisive role and controls the 
amount of stretching of the chains. The free energy per chain (and in units of ksT) is then 



F ~ 



a 



+ V2 




(12) 



The equilibrium height is obtained by minimizing Eq. (12) with respect to /i, and the result 
is 



where the numerical constants have been added for numerical convenience in the following 
considerations. The vertical size of the brush scales linearly with the polymerization index iV, 
a clear signature of the strong stretching of the polymer chains, as was originally obtained by 
Alexander [21]. At the overlap threshold, p* ~ a~'^N~^/^, the height scales as ho ~ iV^/^, and 
thus agrees with the scaling of an unperturbed chain radius in a good solvent, Eq.(5), as it 
should. The simple scaling calculation predicts the brush height h correctly in the asymptotic 
limit of long chains and strong overlap. It has been confirmed by experiments [9, 8, 11] and 
computer simulations [46, 47]. 

3.2 Mean- field calculation 

The above scaling result assumes that all chains are stretched to exactly the same height, 
leading to a step-hke shape of the density profile. Monte-Carlo and numerical mean-field 
calculations confirm the general scaling of the brush height, but exhibit a more rounded 
monomer density profile which goes continuously to zero at the outer perimeter [46]. A big 
step towards a better understanding of stretched polymer systems was made by Semenov [23], 
who recognized the importance of classical paths for such systems. 

The classical polymer path is defined as the path which minimizes the free energy, for 
given start and end positions, and thus corresponds to the most likely path a polymer can 
take. The name follows from the analogy with quantum mechanics, where the classical mo- 
tion of a particle is given by the quantum path with maximal probability. Since for strongly 
stretched polymers the fluctuations around the classical path are weak, it is expected that a 
theory that takes into account only classical paths, is a good approximation in the strong- 
stretching limit. To quantify the stretching of the brush, let us introduce the (dimensionless) 
interaction parameter /3, defined as 



where ho is the brush height according to Alexander's theory, compare Eq. (13). The pa- 
rameter P is proportional to the square of the ratio of the Alexander prediction for the brush 
height, ho, and the unperturbed Gaussian chain radius Ro ~ aN^^'^, and, therefore, is a 
measure of the stretching of the brush based on the scahng prediction. We will later see 
that the actual stretching of the briish is not correctly described by the scaling result for 
short or weakly interacting chains. Constructing a classical theory in the infinite-stretching 
limit, defined as the limit /3 — > cxd, it was shown independently by Milner et al. [48] and 
Skvortsov et al. [49] that the resulting monomer density profile p{z) depends only on the 




(13) 




(14) 
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vertical distance z from the grafting surface and has in fact a parabolic profile. Normalized 
to unity, the density profile is given by [48, 49] 

The brush height, i.e., the value of z for which the monomer density becomes zero, is 
given by z* — {d/n'^Y^^ho and is thus proportional to the scaling prediction for the brush 
height, Eq.(13). The parabolic brush profile has subsequently been confirmed in computer 
simulations [46, 47] and experiments [9] as the limiting density profile in the strong-stretching 
limit, and constitutes one of the cornerstones in this field. Intimately connected with the 
density profile is the distribution of polymer end points, which is non-zero everywhere inside 
the brush (as we will demonstrate later on), in contrast with the original scaling description 
leading to Eq. (13). 

However, deviations from the parabolic profile become progressively important as the 
length of the polymers N or the grafting density pa decreases. In a systematic derivation of 
the mean-field theory for Gaussian brushes [50] it was shown that the mean-field theory is 
characterized by a single parameter, namely the stretching parameter /3. In the limit /3 — > oo, 
the difference between the classical approximation and the mean-field theory vanishes, and 
one obtains the parabolic density profile. For finite f3 the full mean-field theory and the 
classical approximation lead to different results and both show deviations from the parabolic 
profile. 

In Fig. 3 we show the density profiles (normalized to unity) for four different values of /3, 
obtained with the full mean-field theory [50]. In a) the distance from the grafting surface 
is rescaled by the scaling prediction for the brush height, ho, and in b) it is rescaled by 
the unperturbed polymer radius Rq. For comparison, we also show the asymptotic result 
according to Eq. (15) as dashed lines. The self-consistent mean-field equations are solved 
in the continuum limit, where the results depend only on the single parameter /3 and direct 
comparison with other continuum theories becomes possible. As f3 increases, the density 
profiles approach the parabolic profile and already for f3 = 100 the density profile obtained 
within mean-field theory is almost indistinguishable from the parabolic profile denoted by 
a thick dashed line in Fig. 3a. What is interesting to see is that for small values of the 
interaction parameter /3, the numerically determined density profiles exhibit a larger brush 
height than the asymptotic prediction. This has to do with the fact that entropic effects, due 
to steric polymer repulsion from the grafting surface, are not accounted for in the infinite- 
stretching approximation. Experimentally, the achievable f3 values are below /? ~ 50, which 
means that deviations from the asymptotic parabolic profile are important. For moderately 
large values of /3 > 10, the classical approximation (not shown here), derived from the 
mean-field theory by taking into account only one polymer path per cud-point position, 
is still a good approximation, as judged by comparing density profiles obtained from both 
theories [50], except very close to the surface. Unlike mean-field theory, the classical theory 
misses completely the depletion effects at the substrate. Depletion effects at the substrate 
lead to a pronounced density depression close to the grafting surface, as is clearly visible in 
Fig. 3. 

Let us now turn to the thermodynamic behavior of a polymer brush. Using the Alexander 
scaling model, we can calculate the free energy per chain by putting the result for the optimal 
brush height, Eq. (13), into the free-energy expression, Eq. (12), and obtain 

F^(3^N{v2Pa/af\ (16) 
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In Figure 4a we show the rescaled free energy per chain within the full mean-field framework 
(solid line) in comparison with the infinite stretching result (dotted horizontal line) and 
including the leading correction due to the finite entropy of the end-point distribution (broken 
line) [50]. In the infinite-stretching limit, i.e. for — > 0, all curves converge. The free 
energy is not directly measurable, but the lateral pressure can be determined using the 
Langmuir film balance technique. The osmotic surface pressure 11 is related to F, the free 
energy per chain, by 

djFp^A) _ dF _2p^ ( dF/f3\ 

In Figure 4b we show the rescaled osmotic pressure obtained within the mean-field approach 
(solid line). In the infinite-stretching limit, one expects a scaling as 11 ~ p^/^ which is shown 
as a dotted fine. One notes that the more accurate mean-field result gives a pressure which 
is strictly larger than the asymptotic infinite stretching result, in agreement with previous 
calculations and experiments [51, 12, 52, 53, 54]. In the presence of excluded- volume corre- 
lations, i.e., when the chain overlap is rather moderate, the scaling of the brush height h, 
Eq.(13), is still correctly predicted by the Alexander calculation, but the prediction for the 
free energy, Eq.(16), is in error. Including correlations [21], the free energy is predicted to 
scale as F ~ Np^J^. leading to a pressure which scales as 11 ~ p]^/^ in the presence of cor- 
relations. However, all these theoretical predictions do not compare well with experimental 
results for the surface pressure of a compressed brush [12] which has to do with the fact that 
experimentally, chains are rather short so that one essentially is dealing with a crossover 
situation between the mushroom and brush regimes. An alternative theoretical method to 
study tethered chains is the so-called single-chain mean-field method [51], where the sta- 
tistical mechanics of a single chain is treated exactly, and the interactions with the other 
chains are taken into account on a mean-field level. This method is especially useful for short 
chains, where fluctuation effects are important, and for dense systems, where excluded vol- 
ume interactions play a role. The calculated profiles and brush heights agree very well with 
experiments and computer simulations. Moreover, these calculations explain the pressure 
isotherms measured experimentally [12] and in molecular-dynamics simulations [55]. 

A further interesting question concerns the behavior of individual polymer paths. As 
was ahcady discussed for the infinite-stretching theories {j3 oo), polymers paths do end 
at any distance from the surface. In the left part of Figure 5 we show mean-field results for 
the rescaled averaged polymer paths which end at a certain distance from the wall for 
j3 = 1, 10, 100 (from top to bottom). Analyzing the polymer paths which end at a common 
distance from the surface, two rather unexpected features are obtained: i) free polymer ends, 
in general, are stretched; and, ii) the end-points lying close to the substrate are pointing 
towards the surface (such that the polymer path first turns away from the grafting surface 
before moving back towards it). In contrast, end-points lying beyond a certain distance from 
the substrate, point away from the surface (such that the paths move monotonously towards 
the surface). As we will explain shortly below, these two features have been confirmed 
in molecular-dynamics simulations [56]. They are not an artifact of the continuous self- 
consistent theory used in Ref. [50] nor are they due to the neglect of fiuctuations. These 
are interesting results, especially since it has been long assumed that free polymer ends are 
unstretched, based on the assumption that no forces act on free polymer ends. The thick solid 
line shows the unconstrained mean path obtained by averaging over all end-point positions. 
Note that the end-point stretching is small but finite for all finite stretching parameters /3. 
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The right panel exhibits the end-point distributions, which, obviously, are finite over the 
whole range of vertical heights. 

3.3 Molecular Dynamics simulations 

We now present results from Molecular Dynamics simulations in which all the chain monomers 
are coupled to a heat bath. The chains interact via the repulsive portion of a shifted Lennard- 
Jones potential with a Lennard- Jones diameter cr, which corresponds to a good solvent situa- 
tion. For the bond potential between adjacent polymer segments we take a FENE (nonlinear 
bond) potential which gives an average nearest-neighbor monomer-monomer separation of 
typically a ~ 0.97a. In the simulation box with a volume L x L x there are 50 (if not 
stated otherwise) chains each of which consists of N + 1 monomers with varying monomer 
number N — 20, 30, 50, as indicated. The box length L in x— and y— direction was chosen 
to give anchoring densities pa from 0.02 cr"^ to 0.17 cr~^. The first monomer of each chain is 
firmly and randomly attached to the grafting surface at z = 0. The mean square end-to-end 
distance of identical free (i.e. unanchored) chains was for the case = 50 found to be 
Rp = ((ftv — ro)^)f/ee = 12.55(7. For more details consult [56]. 

Figure 6 shows the behavior of the reduced monomer density p{z)Rf/N at increasing 
anchoring density. The stretching of the chains with increasing surface coverage, which is 
due to the repulsion between monomers, is evident. This plot has to be compared with 
Figure 3b, where the same type of rescaling has been used. However, note that at this point, 
direct and quantitative comparison is not possible, since it is a priori not clear which value 
of the interaction parameter (5 in the self-consistent calculation corresponds to which set of 
simulation parameters cr, A^, pa. 

The theoretically predicted scaling law of the brush height Hq ~ N p]J^ has been confirmed 
in several simulations [47]. Provided Pa is above the critical overlap density p* ~ N~^^^, 
the brush height, as measured by the first moment of the monomer density distribution 
{z) — Jzp{z)dz/ fp{z)dz, should approach the predicted scaling form. Figure 7 shows this 
behavior in an appropriate scaling plot of the average brush height {z) and the ^-component 
of the radius of gyration Rg = {J2f=oi'i^i — '^cmf'/{N + 1))^^^ where Vcm is the position of 
the center of mass. For the scaling plot the brush height and radius of gyration arc divided 
by the scaling prediction ~ N p]J^a^/^ and plotted as a function of the analogue of the 
interaction parameter, namely ho/Rp ~ Npl^^a^^^ /Rp, where a and Rp are determined from 
the simulation in the absence of a grafting wall. In this way we obtain an universal crossover 
point N pl/^a^^"^ / Rp ^ 1.4. Thus, for the chain length A^=50, the case we discussed most, 
one can be sure to reach the asymptotic scaling regime for the larger grafting densities. 

A general problem when comparing experimental, simulation and analytical results among 
each other is that the different parameters have to be matched in a meaningful way. One 
such way is based on the relative stretching of chains in a brush. In Figure 8 we plot in a) 
simulation results for the stretching ratio 

7 = {z,{pa))/{Ze{pa = 0)) (18) 

which is the averaged end-point height for a finite grafting density divided by the end- 
point height for vanishing grafting density. The stretching ratio is plotted as a function 
of the scaled anchoring density {N p]J^a^/^ / RpY which is the analogue of the interaction 
parameter (3. It is seen that the different chain lengths and grafting densities scale quite 
nicely. Clearly, for short chains or low grafting densities, that is in the mushroom regime, the 
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stretching ratio 7 approaches unity. The highest stretching ratio reached in the simulations 
is 7 2.7. In Figure 8b we show mean- field results for the stretching ratio 7 as a function 
of the interaction parameter /3. The general shape of the curve is similar. Matching the 
stretching ratios of the simulation and the mean-field calculation allows to determine the 
effective (3 value corresponding to a particular simulation run. As an example, for the most 
highly stretched simulation with 7 fti 2.7, we estimate an effective interaction parameter of 
/3 25. A similar matching procedure is possible with experimental data. 

In agreement with the mean- field profiles shown in Figure 3, the monomer density ob- 
tained from the simulations in Figure 6 decays smoothly to zero at distances far from the 
anchoring plane. However, comparing the shape of the profiles in more detail, one notes a re- 
gion over which the simulation profiles become rather fiat for anchoring densities pa > O.lcr^^, 
in contrast to the mean-field results. According to Figure 8, this corresponds to stretching 
factors 7 > 2 or /3 > 10. This discrepancy can be traced back to the neglect of higher-order 
virial terms in the analytical theory, which beome important at elevated monomer densities. 

Wc also looked at individual polymer paths within the simulation. The scaled average 
polymer paths with a given fixed end-point distance Zg, {zz^{n)) / (ze) , are shown in Figure 9 
for three typical anchoring densities corresponding to stretching factors 7 ~ 1.3, 1.8, and 
2.7, from top to bottom. To get an idea of their relative weight we plot the polymer paths 
together with the normalized density of free ends Pe{z) (shown in the right panel). The 
trajectories are very similar to those predicted by the mean- field theory (see Figure 5). 
Paths which end far from the anchoring surface are stretched through their entire length, 
including the free end point. The paths which end at the outer rim of the distribution are 
almost uniformly stretched and appear almost as a straight line. On the other hand, paths 
which end close to the anchoring surface arc non-monotonic, first moving away from the wall, 
reaching a maximal distance, and then turning back towards the anchoring surface. Except 
at the maximum, all paths are stretched everywhere, including the end point. The straight 
broken lines are added for comparison and denote the maximally stretched possible path, 
which consists of a fully oriented and straight polymer. As can be seen from Figure 5 as well 
as from Figure 9, the end-point stretching, proportional to (dz/ds) or {dz/dn), respectively, 
is positive for some particular paths and negative for others, so that the average, plotted as 
thick line in the Figures, is typically quite small. This reconciles the present results with the 
infinite-stretching results by Milner, Witten and Gates [48] and Zhulina and co-workers [49] , 
where the chains were assumed to be unstretched at their ends regardless of their end-point 
position: this assumption turns out to be acceptable on average if the stretching of the chains 
is large, i.e., for large 

3.4 Additional effects 

As we described earlier, the main interest in end-adsorbed or grafted polymer layers stems 
from their ability to stabilize surfaces against van-der-Waals attraction. The force between 
colloids with grafted polymers is repulsive if the polymers do not adsorb on the grafting 
substrates [57], that is, in the absence of polymer bridges and loops. A stringent test 
of brush theories was possible with accurate experimental measurements of the repulsive 
interaction between two opposing grafted polymer layers using a surface force apparatus [8]. 
The resultant force could be fitted very nicely by the infinite-stretching theory of Milner et 
al. [58]. It was also shown that polydispersity effects, as appear in experiments, have to be 
taken into account theoretically in order to obtain a good fit of the data [59]. 
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So far we assumed that the polymer grafted layer is in contact with a good solvent. 
In this case, the grafted polymers try to minimize their mutual contacts by stretching out 
into the solvent. If the solvent is bad, the monomers try to avoid the solvent by forming a 
collapsed brush, the height of which is considerably reduced with respect to the good-solvent 
case. It turns out that the collapse transition, which leads to phase separation in the bulk, 
is smeared out for the grafted layer and does not correspond to a true phase transition [60] . 
The height of the collapsed layer scales linearly in PaN , which reflects the constant density 
within the brush, in agreement with experiments [61]. Some interesting effects have been 
described theoretically [62] and experimentally [61] for brushes in mixtures of good and bad 
solvent, which can be rationalized in terms of a partial solvent demixing. 

For a thcta solvent {v2 = 0) the relevant interaction is described by the third-virial 
coefficient; using a simple Alexander approach similar to the one leading to Eq. (13), the 
brush height is predicted to vary with the grafting density as /i ~ p^^, in agreement with 
computer simulations [63]. 

Up to now we discussed planar grafting layers. It is of much interest to consider the case 
where polymers are grafted to curved surfaces. The first study taking into account curvature 
effects of stretched and tethered polymers was done in the context of star polymers [24]. 
It was found that chain tethering in the spherical geometry leads to a universal density 
profile, showing a densely packed core, an intermediate region where correlation effects are 
negligible and the density decays as p(r) ~ 1/r, and an outside region where correlations are 
important and the density decays as p ~ r~^/^, with the radial distance denoted by r. These 
considerations were extended using the infinite-stretching theory of Milner et al. [64], self- 
consistent mean-field theories [65], and molecular-dynamics simulations [66]. Of particular 
interest is the behavior of the bending rigidity of a polymer brush, which can be calculated 
from the free energy of a cylindrical and a spherical brush and forms a conceptually simple 
model for the bending rigidity of a lipid bilayer [67]. 

The behavior of a polymer brush in contact with a polymeric solvent, consisting of chem- 
ically identical but somewhat shorter chains than the brush, had been first considered by de 
Gennes [22]. A complete scaling description has been given only recently [68]. One distin- 
guishes different regimes where the polymer solvent is expelled to various degrees from the 
brush. A somewhat related question concerns the behavior of two opposing brushes brought 
closely together, and separated by a solvent consisting of a polymer solution [69, 70]. Here 
one distinguishes a regime where the polymer solution leads to a strong attraction between 
the surfaces via the ordinary depletion interaction, but also a high polymer concentration 
regime where the attraction is not strong enough to induce colloidal fiocculation. This 
phenomenon is called colloidal rcstabilization [69]. 

Considering a mixed brush made of mutually incompatible grafted chains, a novel transi- 
tion to a brush characterized by a lateral composition modulation was found [71] . Even more 
complicated spatial structures are obtained with grafted diblock copolymers [72]. Finally, 
we would like to mention in passing that these static brush phenomena have interesting 
consequences for the dynamic properties of polymer brushes [73] . 

4 Charged grafted polymers 

Brushes can also be formed by charged polymers which are densely end-grafted to a sur- 
face; the resulting charged brush shares many of the features discussed above for neutral 
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brushes, but qualitatively new properties emerge due to the presence of charged monomers 
and counterions in the brush. Charged brushes have been the focus of numerous theoretical 
[74]-[82] and experimental [83]-[87] studies. They serve as an efficient mean for preventing 
colloids in polar media (such as aqueous solutions) from flocculating and precipitating out 
of solution. This stabilization arises from steric (entropic) as well as electrostatic repulsion. 
A strongly charged brush is able to trap its own counter-ions and generates a layer of locally 
enhanced salt concentration [76] . It is thus less sensitive to the salinity of the surrounding 
aqueous medium than a stabilization mechanism based on pure electrostatics (i.e. without 
polymers). Compared to the experimental knowledge about uncharged polymer brushes, 
less is known about the scaling behavior of PE brushes. The thickness of the brush layer 
has been inferred from neutron-scattering experiments on end-grafted polymers [83] and 
charged diblock-copolymers at the air- water interface [85]. 

Theoretical work on PE brushes was initiated by the works of Miklavic and Marcelja 
[74] and Misra et al. [75]. In 1991, Pincus [76] and Borisov, Birshtein and Zhulina [77] 
presented scaling theories for charged brushes in the so-called osmotic regime, where the 
brush height results from the balance between the chain elasticity (which tends to decrease 
the brush height) and the repulsive osmotic counter-ion pressure (which tends to increase the 
brush height). In later studies, these works have been generalized to poor solvents [78] and 
to the regime where excluded volume effects become important, the so-called quasi-neutral 
or Alexander regime [81]. 

4.1 Scaling approach 

An analytical theory for polyelectrolyte brushes relies on a number of simplifying assump- 
tions. The full theoretical problem is intractable because the degrees of freedom of the 
polymer chains and the counterions are coupled by steric and long-ranged Coulomb inter- 
actions. It is important to note that charged polymers by itself are not fully understood, 
therefore quite drastic simplification are needed to tackle the more complicated system of 
polyelectrolytes end-grafted to a surface. Firstly, we will concentrate on polymeric systems 
with counterions and only briefly mention the effects of added salt towards the end of this 
section, which has been discussed extensively in the original literature. Secondly, we will 
write the total free energy per unit area and in units of ksT, 

F = Fp,, + Fi,^ + Fi^t (19) 

as a sum of separate contributions from the polyelectrolytes. Fpoi, contributions from the 
counterions, F[o^, and an electrostatic interaction term Fi^t which couples polymers and 
counterions. The schematic geometry of the brush system is visualized in Figure 10: we 
assume that the charged brush is characterized by two length scales. The polymer chains 
are assumed to extend to a distance h from the grafting surface, the counterions in general 
form a layer with a thickness of d. Two different scenarios emerge: The counter-ions can 
either extend outside the brush, d ^ h, as shown in Figure 10a, or be conflned inside the 
brush, d ^ h as shown in Figure 10b. As we demonstrate now, case b) is indicative of 
strongly charged brushes, and thus applicable to most experiments done on charged brushes 
so far, while case a) is typical for weakly charged brushes. 

We recall that the grafting density of PE's is denoted hy pa, q is the counter- ion valency, 
the polymerization index of grafted chains, and / the charge fraction. The counterion free 
energy contains entropic contributions (due to the confinement of the counterions inside a 
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layer of thickness d) and also energetic contributions which come from interactions between 
counterions. In previous theories, a low-density expansion for the interaction part was used. 
We remind that the second-virial interaction is important for neutral brushes and is the 
driving repulsive force balancing the chain elasticity. For charged brushes, on the other 
hand, the leading term of the electrostatic correlation energy (which shows fractional scaling 
with respect to the charge density) is attractive and has been shown to lead to a collapse 
of the polyelcctrolyte brush for large Bjerrum lengths[88, 89]. We will not pursue this 
collapse transition further in this review paper but solely concentrate on entropic and steric 
counterion effects, which is an acceptable approximation for mono-valent counterions.. In 
the present analysis we use a free-volume approximation very much in the spirit of the 
van-der-Waals equation for the liquid-gas transition. We include the effective hard-core 
volume of a single polyelcctrolyte chain, which we call V, and which reduces the free volume 
that is available for the counterions. This free volume theory therefore takes the hard-core 
interactions between the polymer monomers and the counterions into account in a non- linear 
fashion and is valid even at large densities in the limit of close-packing. Compared to that, 
the excluded-volume interaction between counterions is small since the monomers are more 
bulky than the counterions and therefore it is neglected. The non-linear entropic free energy 
contribution of the counterions reads per unit area 



In the limit of vanishing polymer excluded-volume, — > 0, one recovers the standard ideal 
entropy expression. As the volume available for the counterions in the brush, which per 
polymer is just d/pa, approaches the self volume of the polymers, V, the free energy expres- 
sion Eq. (20) diverges, that means, the entropic prize for that scenario becomes infinitely 
large. The excluded volume of the polymers is roughly independent of the polymer brush 
height, and can be written in terms of the efi^ective monomer hardcore diameter a^s and the 
polymer contour length aN as V ~ aNa'^g. This leads to the final expression 



The polymer free energy Fpoi per unit area can in principle be derived from the free en- 
ergy per chain used for neutral brushes. Equation (12), by multiplying with the grafting 
density. It contains two contributions, one due to the chain elasticity and the other due 
to monomer-monomer interactions. Note that a logarithmic entropy term as for the coun- 
terions is missing, since the translational degrees of freedom of the monomers are lost due 
to chain connectivity. We recall that the entropy of chain stretching was accounted for by 
linear elasticity theory, with a spring constant proportional to 1/Rq in units of the thermal 
energy. This linear expression for the chain stretching is reliable for the neutral case, where 
indeed chain stretching is typically quite mild. For charged chains, the stretching is much 
stronger and we have to consider a more detailed model for the chain elasticity. For a freely 
jointed chain, which is a good model for synthetic polymer chains, the entropy loss due to 
stretching can be calculated exactly [90]. We only need here the asymptotic expressions for 
weak and for strong stretching, which read (per unit area) 




(20) 
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and are proportional to the grafting density pa- The contour length of the fully stretched 
chain is Na. The weak-stretching term is the standard term used in previous scahng mod- 
els. For the highly stretched situations encountered in fully charged brushes, the strong- 
stretching term is typically more appropriate and leads to a few changes in the results as 
will be explained further below. The energetic contribution to the polymer free energy can 
be expressed as the second virial contribution, arising from steric repulsion between the 
monomers (contributions due to counter ions are neglected). Throughout this section, the 
polymers are assumed to be in a good solvent (positive second virial coefficient > 0) . The 
contribution thus reads 

^' - 5""^ (^)" • <^^' 

An additional electrostatic component to the polymer interaction term is typically unim- 
portant since the counterions strongly screen any Coulomb interactions [90]. Finally, an 
electrostatic interaction between polymers and counterions Fjnt occurs if the PE brush is not 
locally electro-neutral throughout the system, as for example is depicted in Fig. 10a. This 
energy is given by 

2t^I^{N f p^f {d - hf 

Fint = ^ ^ . (24) 

This situation arises in the limit of low charge, when the counter-ion density profile extends 
beyond the brush layer, i.e. d > h. 

The different free energy contributions lead, upon minimization with respect to the two 
length scales h and d, to different behaviors. Let us first consider the weak charging limit, 
i.e. the situation where the counter- ions leave the brush, d > h. In this case, minimization 
of Fion + -^int with respect to the counter-ion height d in the limiting case of vanishing brush 
height {h — 0) and monomer volume {a — 0) leads to 

d = ^ ,1,, = 3Agc (25) 

which has the same scaling as the so-called Gouy-Chapman length Xqc- This length scale 
is a measure for the average height of the diffuse layer of q'-valent counter-ions adsorbed 
at a surface with effective surface charge density Nfpa and has been determined within 
simulations and ficld-thcory[91] but can also be obtained with more coarse-grained scaling 
arguments, as demonstrated here. Balancing now the polymer stretching energy Fp^j and 
the electrostatic energy Fjnt one obtains the so-called Pincus brush height 

h^N^Paa^Bf , (26) 

which results from the electrostatic attraction between the counter-ions and the charged 
monomers. One notes the peculiar dependence on the polymerization index A'^. In the 
limit oi d !^ h, the PE brush can be considered as neutral and the electrostatic energy 
vanishes. There are two ways of balancing the remaining free energy contributions. The first 
is obtained by comparing the osmotic energy of counter-ion confinement, -Fion, in the limit 
when d = h and for vanishing polymer volume, with the polymer stretching term, F^^i, in 
the weak stretching limit, leading to the height 
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constituting the linear osmotic brush regime. The main assumption here is that all counte- 
rions stay strictly localized inside the brush, which will be tested later by comparison with 
computer simulations. 

Finally comparing the second- virial term for the counterion interactions, -Fpoi, with the 
polymer stretching energy in the weak-stretching limit, F^^^, one obtains the same scaling 
behavior as the neutral brush [21, 22], compare Eq.(13). Comparing the brush heights in all 
three regimes we arrive at the phase diagram shown in Fig. 11. The three scaling regimes 
coexist at the characteristic charge fraction 



qv2 



1/3 



and the characteristic grafting density 



4.2 Non- linear effects 

The scaling relations for the brush height and the crossover boundaries between the various 
regimes constitute the simplest approach towards charged brushes. We pointed out already 
a few limitations of the presented results, which have to do with non-linear stretching and 
finite- volume effects. 

Computer simulations provide an excellent mean to test those scaling relations. Extensive 
molecular dynamics simulations have been performed recently for polyelectrolyte brushes at 
various grafting densities and charge fractions, both at strong and intermediate electrostatic 
couplings[92, 93]. In these simulations, a freely-jointed bead-chain model is adopted for 
charged polymers end-grafted onto a rigid surface. The counterions are explicitly modeled 
as charged particles where both counterions and charged monomers are univalent and inter- 
act with the bare Coulomb potential Eq. (8). The strength of the Coulomb interaction is 
controlled by the Bjerrum length is- No additional electrolyte is added. The short-range 
repulsion between all particles is modeled by a shifted Lennard- Jones (LJ) potential, charac- 
terized by the hard-core diameter a, as was done for the neutral brush simulations in Section 
3.3. The simulation box is periodic in lateral directions and finite in the z-direction normal 
to the anchoring surface at z—0. We apply the techniques introduced by Lekner and Sperb to 
account for the long-range nature of the Coulomb interactions in a laterally periodic system 
[92] . To study the system in equilibrium we use stochastic molecular simulation at constant 
temperature. 

Simulated density profiles of monomers and counterions of the system in normal direc- 
tion are shown in Figure 12 for the fully charged brush at several grafting densities and for 
a Coulomb coupling characterized by is/o' = 1, which is close to the coupling of mono- 
valent ionic groups in water. As seen, both monomers and counterions follow very similar 
nearly-step-like profiles with uniform densities inside the brush, which increases with grafting 
density. These data show that the counterions are mostly confined in the brush layer and 
that the electroneutrality condition is satisfied locally. The simple explanation is that the 
Gouy-Chapman length is indeed very small and of the order or smaller than the monomer 
diameter. One can observe that the polyelectrolyte chains are stretched up to about 70 % 
of their contour length, which is roughly L ^ Na = 30a, and thus their elastic behavior is 
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far beyond the linear regime. Therefore, within the chosen range of parameters, the sim- 
ulated brush is in the strong-charging (i.e., all counterions are confined within the brush) 
and strong-stretching limits; as we will show below, it exhibits a non-linear osmotic scaling 
behavior, so the linear scaling description of the last section, leading to Eq. (27), has to be 
slightly modified. The average height of end-points of the chains, (Zg), is one of the quantities 
which can be directly measured in the simulations and is shown in Figure 13 together with 
the predictions of the hnear scaling theory, Eq. (27), denoted by (a), and the non-hnear 
scaling predictions, to be explained further below. It is observed that the simulated brush 
height (solid circles) varies slowly with the grafting density, contrary to the prediction of the 
linear scaling theory, Eq.(27), but in agreement with recent experimental results[94, 95]. 

We now present a scaling theory that incorporates non-linear elastic and osmotic ef- 
fects. In order to bring out the physics most clearly, the derivation will be done in two 
consecutive steps. First, in the strongly-stretched osmotic brush regime, one chooses the 
strong-stretching version of the chain-stretching entropy in Equation (22) and balances it 
with the counterion entropy. Equation (21), assuming vanishing polymer self volume, o" = 0, 
and assuming equal heights for the brush and counterion layers, d = h. The result is 

Naqf 

which is the large-stretching analogue of Equation (27). The maximal stretching predicted 
from this equation is obtained for f=l and corresponds to a vertical chain extension corre- 
sponding to 50 % of the contour length. For comparison, both expressions (27) and (30) are 
shown in Figure 13 as dashed lines (a) and (b) respectively. Still, the overall stretching is 
considerably smaller than what is observed in experiments and simulations, and it transpires 
that something is missing in the above scaling description. It has to do with the entropic 
pressure which increases as the volume within the brush is progressively more filled up by 
monomers and counterions. Note that the brush height in Equations (27) and (30) does 
not depend on the grafting density, in vivid contrast to the simulation results displayed in 
Figure 13. 

In the non-linear osmotic brush regime we combine the high-stretching (non-linear) ver- 
sion of the chain elasticity in Equation (22) with the non-linear entropic effects of the counte- 
rions due to the finite volume of the polymer chains, i.e. we choose a finite polymer diameter 
a in Equation (21). The final result for the equilibrium brush height is 

Na{qf + a^spa) 
' OT^' ^''^ 

which in the limit of maximal grafting density, that is close packing, pa l/'^eff; reaches the 
maximal value stretching value, h — > Na, as one would expect: Compressing the brush lat- 
erally increases the vertical height and finally leads to a totally extended chain structure. In 
Figure 13, we compare expression Eq. (31) with the simulation results for the brush height as 
a function of grafting density. Note that we have used a^^ = 2a^. This choice corresponds to 
an approximate two-dimensional square-lattice packing of monomers and counterions on two 
interpenetrating sublattices. As can be seen, the scaling prediction, Eq. (31), qualitatively 
captures the slow brush height increase with grafting density, as observed in the simulations. 
The deviations between simulation data and the non-linear osmotic brush prediction are un- 
problematic since the effective hard-core diameter of monomers in the simulation is difficult 
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to determine precisely and might as well be treated as a fitting parameter. Deviations might 
also be explained by considering higher-order effects, such as lateral inhomogeneity of the 
counterion distribution around the brush chains and intermediate-stretching elasticity of the 
chains, that go beyond the present scahng analysis and have been considered in Ref.[90]. 

4.3 Additional effects 

For large values of the charge fraction / and the grafting density Pa, and for large Coulomb 
coupling, it has been found numerically that the brush height does not follow any of the 
scaling laws discussed here [92]. This has been recently rationalized in terms of another 
scaling regime, the collapsed regime. In this regime one finds that correlation and fluctuation 
effects, which are neglected in the discussion in this section, lead to a net attraction between 
charged monomers and counter- ions [88, 89]. 

If salt is present in the solution, counterions as well as co-ions do penetrate into the brush, 
which leads to additional screening of the Coulomb repulsion inside the brush. The amount 
of this screening, and the stretching of the polyectrolyte chains, are now also controlled 
by the bulk salt concentration. Since the additional salt screening weakens the sweUing 
of the brush caused by the counterion osmotic pressure, salt leads to a brush contraction 
for sufficiently high salt concentration according to ~ CgJ/^[76, 77]. The threshold salt 
concentration above which the brush contraction sets in is given by the salt concentration 
which equals the counterion concentration inside the brush. This means that the higher the 
grafting density (and consequently the higher the internal counterion concentration in the 
osmotic brush regime), the larger the salt concentration neceassary to see any salt effects at 
all. 

Another way of creating a charged brush is to dissolve a diblock copolymer consisting of 
a hydrophobic and a charged block in water. The diblocks associate to form a hydrophobic 
core, thereby minimizing the unfavorable interaction with water, while the charged blocks 
form a highly charged corona or brush [96]. The micelle morphology depends on different 
parameters. Most importantly, it can be shown that salt acts as a morphology switch, 
giving rise to the sequence spherical, cylindrical, to planar micellar morphology as the salt 
concentration is increased [96]. Theoretically, this can be explained by the entropy cost of 
counter-ion confinement in the charged corona [97]. The charged corona can be studied by 
neutron scattering [98] or atomic-force microscopy [99] and gives information on the behavior 
of highly curved charged brushes. 

5 Concluding Remarks 

In this short review about neutral and charged polymers that arc terminally grafted with one 
end to a surface we tried to explain and contrast the most important theoretical methods used 
for their understanding and description. The theories used for brushes are quite different 
from ordinary polymer modelling since the statistics of the grafted layer depends crucially 
on the fact that the chain is not attracted to the surface but is forced to be in contact 
to the surface since one of its ends is chemically or physically bonded to the surface. We 
review scaling concepts, mean-field theory and simulation techniques that give information 
on brushes on different levels of detail and accuracy. Scahng concepts yield the qualitative 
features of brushes, i.e. the brush height and its free energy as a function of the main 
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parameters and without reliable prefactors. Mean-field or self- consistent theories allow to 
construct density distributions on a coarse grained level, but correlation effects are completely 
missed. Simulation techniques give in principle exact numerical results for the case of so- 
called primitive models where the solvent is neglected and monomers and ions are modeled 
as charged soft spheres. Current open questions concern the structure of water in such dense 
systems: is water still described by a continuous medium with the bulk dielectric constant? 
Secondly, all approaches neglect the varying polarizability of the monomers and ions, which 
again give rise to pronounced image-charge effects in the brush geometry. Lastly, most brush 
systems have not reached their true equilibrium structure, so that we are in essence dealing 
with an non-equilibrium system. How do we describe such non-equilibrium systems, and 
what governs the approach towards equilibrium including hydrodynamic and local friction 
effects? We believe that the triad of scaling, field-theoretic, and simulation techniques will 
also allow to gain understanding of these more complicated issues in the near future. 
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a) 




Figure 1: For grafted chains, one distinguishes between: a) the mushroom regime, where 
the distance between chains, p^^^^, is larger than the size i? of a polymer coil; and, b) the 
brush regime, where the distance between chains is smaller than the unperturbed coil size. 
Here, the chains are stretched away from the surface due to repulsive interactions between 
monomers. The brush height h scales linearly with the polymerization index, h ~ N. and 
thus is larger than the unperturbed coil radius R which in the good-solvent regime scales 
according to Flory as i?^ ~ aN'^. 
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Figure 2: Top: Freely-jointed chain (FJC) model, where N bonds of length a are connected 
to form a flexible chain with a certain end-to-end distance R. Bottom: In the simplified 
model, appropriate for more advanced theoretical calculations, a continuous line is governed 
by some bending rigidity or line tension. This continuous model can be used when the 
relevant length scales are much larger than the monomer size. 
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Figure 3: Self-consistent mean-field results for the density profile (normalized to unity) of a 
brush for different values of the interaction parameter /?. In a) the distance from the grafting 
surface is rescaled by the scaling prediction for the brush height, ho, and in b) it is rescaled 
by the unperturbed polymer radius Rq. As /3 increases, the density profiles approach the 
parabolic profile (shown as dashed lines). 



27 




0.2 0.4 0.6 0.8 1 '0 0.2 0.4 0.6 0.8 1 



1/P 1/P 

Figure 4: a) Mean-field result (solid line) for the rescaled brush free energy per polymer as a 
function of the inverse interaction parameter The infinite stretching result is indicated 
by a horizontal dotted line, the broken straight line denotes the infinite stretching result 
with the leading correction due to the finite end-point distribution entropy, b) Rcscalcd 
lateral pressure within mean-field theory (solid line) compared with the asymptotic infinite- 
stretching result (dotted line). 
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Figure 5: Left: Mean- field results for the rescaled averaged polymer paths which end at a 
certain distance Zf. from the wall for f3 = 1, 10, 100 (from top to bottom), corresponding 
to stretching values of 7 = 1.1, 1.9,5.6. The thick solid line shows the unconstrained mean 
path obtained by averaging over all end-point positions. Note that the end-point stretching 
is small but finite for all finite stretching parameters /3. Right: End-point distributions. 
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Figure 6: Simulation results for the normalized monomer number density, p{z)Rf / N pa, as a 
function of the scaled distance from the grafting surface z/Rp for anchored chains of length 
= 50 and grafting densities paCr"^ = 0.02, 0.04, 0.06, 0.09, and 0.17 (from top to bottom). 
Note that Rp is determined within the simulation for a single, free polymer chain. 
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Figure 7: Simulation results for the rescaled average monomer height h/{Npl^^a^/^) (top) 
and the vertical component of the radius of gyration Rgz/ {N p]J^a^/^) (bottom) versus the 
scahng variable N{paa^Y/^ /Rp for different chains lengths as indicated in the figure. The 
open/filled symbols correspond to extensible and non-extensible chains, respectively, which 
is controlled within the simulation by the strength of the FENE bond potentials. As one 
can see, for large values of the parameter combination N{paa^Y^^ /Rp the data saturate at 
a plateau and are thus in agreement with the scaling prediction. 
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Figure 8: Effective stretching factor 7 = (ze(pa)/-2e(0)) as obtained a) witliin tlie simulation 
as a function of tlie scaled anchoring density {Npl/^a^^^ / RpY and b) from the self-consistent 
field theory as a function of the interaction parameter /3. The comparison between the 
simulation results and the analytical results allows to determine the effective /3 parameter 
of a particular simulation set. 
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Figure 9: Left: Simulation results for the rescaled averaged polymer paths which end at a cer- 
tain distance Ze from the wall for pacr^ = 0.02, 0.06, 0.17 (from top to bottom), corresponding 
to stretching values of 7 = 1.3, 1.8, 2.7. The thick solid line shows the unconstrained mean 
path obtained by averaging over all end-point positions. Right: End-point distributions. All 
data are obtained for simulations with 50 chains consisting of A?" = 50 monomers each. 
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Figure 10: Schematic PE brush structure. In a) we show the weakly- charged hmit where the 
counter-ion cloud has a thickness d larger than the thickness of the brush layer, h. In b) we 
show the opposite case of the strongly-charged limit, where all counter- ions are contained 
inside the brush and a single length scale d ^ h exists. 
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Figure 11: Scaling diagram for PE brushes on a log-log plot as a function of the grafting 
density pa and the fraction of charged monomers /. Featured are the Pincus-brush regime, 
where the counter-ion layer thickness is much larger than the brush thickness, the osmotic- 
brush regime, where all counter- ions are inside the brush and the brush height is determined 
by an equilibrium between the counter-ion osmotic pressure and the PE stretching energy, 
and the neutral-brush regime, where charge effects are not important and the brush height 
results from a balance of PE stretching energy and second-virial repulsion. The power law 
exponents of the various lines are denoted by numbers. 
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Figure 12: Density profiles of monomers (open symbols) and counterions (filled symbols) as 
a function of the distance from the anchoring surface. Shown are profiles for fully charged 
brushes of 36 chains ol N — monomers at grafting densities (from bottom to top) paC^ = 
0.042 (circles), 0.063 (squares), 0.094 (diamonds), and 0.12 (triangles). As can be clearly 
seen, the counterion stay inside the brush for all considered grafting densities and the local 
electroneutrality condition is satisfied very well. 
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Figure 13: The circles show simulation data for the brush height as a function of grafting 
density for polyelectrolyte chains oi N — 30 monomers (contour length L — 30a) with charge 
fraction of / = 1. The Bjerrum length is is — cr, which corresponds to an intermediate 
coupling strength quite relevant for fully charged chains in water, and all ions are univalent. 
The dotted lines (a) and (b) show the linear scaling predictions, Eqs. (27) and (30), with 
Gaussian and non-linear elasticity respectively. The solid line shows the non-linear scaling 
prediction according to Eq. (31), which includes non-hnear elasticity as well as the finite 
excluded volume of monomers and counterions. 
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